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ABSTRACT 

We derive radial mass profiles of four strong lensing selected clusters which 
show prominent giant arcs (Abell 1703, SDSS J1446+3032, SDSS J1531+3414, 
and SDSS J2111— 0115), by combining detailed strong lens modeling with weak 
lensing shear measured from deep Subaru Suprime-cam images. Weak lensing sig- 
nals are detected at high significance for all four clusters, whose redshifts range 
from z = 0.28 to 0.64. We demonstrate that adding strong lensing information 
with known arc redshifts significantly improves constraints on the mass density 
profile, compared to those obtained from weak lensing alone. While the mass 
profiles are well fitted by the universal form predicted in iV-body simulations of 
the A-dominated cold dark matter model, all four clusters appear to be slightly 
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more centrally concentrated (the concentration parameters c v - 1T ~ 8) than theoret- 
ical predictions, even after accounting for the bias toward higher concentrations 
inherent in lensing selected samples. Our results are consistent with previous 
studies which similarly detected a concentration excess, and increases the total 
number of clusters studied with the combined strong and weak lensing technique 
to ten. Combining our sample with previous work, we find that clusters with 
larger Einstein radii are more anomalously concentrated. We also present a de- 
tailed model of the lensing cluster Abell 1703 with constraints from multiple 
image families, and find the dark matter inner density profile to be cuspy with 
the slope consistent with —1, in agreement with expectations. 

Subject headings: dark matter — galaxies: clusters: individual (Abell 1703, SDSS 
J1446+3032, SDSS J1531+3414, SDSS J2111— 0115) — gravitational lensing 



Introduction 



The current standard model of structure formation is successful i n explaining variou s 
cosmological observations, such as the cosmic microwave background (jSpergel et al.ll2003l ). 
the clustering pattern of galaxies (lEisenstein et al.ll2005l ). the Lya forest fluctuations in 



the intergalactic med i um ( McDonald et al. 



20051 ) . and the abundance of clusters of galaxies 



( IGladders et al.ll2007l ; IVikhlinin et al.ll2009l ). In the theory, the growth of structure is driven 
by the gravitational instability of dark matter. The initial density fluctuations are nearly 
scale -invariant and Gaussi an, as inferred from the cosmic microwave background observations 



: e.g.. lKomatsu et al.ll2009l ). The spectrum of the fluctuation, together with the "cold" ansatz 
of dark matter, suggests that astronomical objects form through the bottom-up assembly 
process. 

Clusters of galaxies serve as one of the most important tests of the standard structure 
formation model. Clusters are the largest virialized objects in the universe. Their high virial 
temperatures suggest that most of baryon remain hot in massive clusters and dissipative 
cooling of baryons is less efficient. Therefore, the density profile of clusters can well be 
approximated by the distribution of dark matter which is sensitive to both the nature of 
dark matter and the dark matter assemble history. In the standard collisionless cold dark 
matter (CDM) scenario, the density profile of dark matter found in iV-bo dy simulations 



has a universal form with its slope progressively shallower toward the center (jNavarro et al. 



19971 . hereafter NF W). Cluster-scale halos are predicted to have rather low-concentration 



mass profiles (e.g., iBullock et al.ll200ll ; iNeto et ajj 120071 ). These are important theoretical 



predictions that should be confronted with observations. 
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Gravitational lensing is suited for this purpose, because it probes the distribution of 
matter directly, regardless of the light distribution. Weak gravitational lensing, which 
takes advantage of small distortions of background galaxies to reconstruct the mass dis- 
tribution, is particularly powerful in extracting outer mass profiles of massive clusters. On 
the other hand, strong lensing provides robust measurements of cluster masses near the 
center. Therefore, it is essential to combine stron g and weak lensing analysis in order to con 



strain densit y profiles over a wide range of radii (Nataraian et al.lll 998. 



2006. 


2008all 




Diego et al. 


Merten et al. 


20091). 



2007 



Hicks et al. 



2002 



2007 



Bradac et al. 



Deb et al 



2008 



Cluster mass profiles grow i n importance, particularly given possible high concentrations 
found in some lens-rich clusters. iBroadhurst et al.l (j2005al ) argued that the radial mass profile 
of Abell 1689 (A1689) from combined strong and weak lensing analysis appears to be sig- 
nificantly m ore centrally concentrated than expected from iV-body simulations of dark mat 



ter (see also Clowe fe Schneider! 200 ll; 



2007; Limousin et al 



2007 



King et al. 12002 ; lHalkola et al.l 12Q06J ; iBardeau et al. 



Umetsu fe Broadhurstll2008l ). S uch high concentrat ions have also 



been suggested in o ther lensing clusters, such as CL0024 (jKneib et al.ll2003l ) and MS2137 
(IGavazzi et al.ll2003l ). A possible resolution i s to consider a triaxial cluster with i ts major 



axis aligned with the line-o f-sight direction (I Clowe et al.l 120041 ; lOguri et al.l 120051 ; iGavazzi 



20051 ; ICorless fe King 120071 ) . Because of t he lensing bias that lens-rich clusters are more 
likely to be oriented with the major axis (IHennawi et al.l 120071 ) . the high conc entration of 



A1689 can in practice be marginally reconciled with the ACDM predi ctions (lOguri et al. 



20051 : IHennawi et al.ll2007l : lOguri fe Blandfordll2009l ; ICorless et al.l 120091 ) . The broad range 
of apparent concentrations expected in the ACDM model, both from the intrinsic scatter of 
concentrations and from additional scatter due to the projection effect, indicates the need 
for t he statistical study of c oncentrations using a well-controlled sample of lensing clusters 



(.e.g. 



Broadhurst et al.ll2008h . 



In this paper, we present new detailed lensing analysis of four clu sters. All the clusters 
are new strong lensing clusters discovered by our survey of giant arcs (IHennawi et al.ll2008l ) 
using a catalog of massive clusters co nstructed from th e imaging and spectroscopic data of 
the Sloan Digital Sky Survey (SDSS; lYork et al.llioooh. We exploit wide field-of-view and 



excellent image quality of the Subaru Suprime-cam (IMiyazaki et al.l 120021 ) to determine the 
outer mass profiles of the clusters with weak lensing technique. The weak lensing signals are 
combined with strong lensing constraints with redshifts of arcs obtained from the Gemini 
telescope. The resulting radial mass profiles of these clusters constrained from strong and 
weak lensing are compared with the ACDM prediction. 



The structure of this paper is as follows. We summarize our cluster sample for analysis 
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and follow-up data in §2j Then we perform strong and weak lensing analyses in §[3] and 
§U respectively. These results are combined in ||5l Our results are summarized in $61 In 
Appendix [AJ we discuss the inner density profile of one of our target clusters, Abell 1703, 
constrained from strong lens modeling. Throughout the paper, we adopt a cosmological 
model with the matter density Qm = 0.26, the cosmological constant = 0.74, and the 
dimensionless Hubble constant h = 0.72. 



2. Cluster Sample and Follow-up Data 



We selected our targets for detailed lensing analysis f rom our sample of le nsing clusters 
with definitive arcs in the SDSS Giant Arc Survey (SGAS; lHennawi et al.ll2008l ). The cluster 
lens sample has been constructed as follows. We take advantage of a sample of the richest 
clusters fro m the 8000 deg 2 photometric data of the SDSS using a red-sequence cluster finding 



algorithm (IGladders fc Yeell2000l . 120051 ) . We then obtained follow-up images of the richest 



clusters with the Wisconsin Indiana Yale NO AO (WIYN) 3.5m telescope, the University of 
Hawaii 88-inch telescope (UH88), and the Nordic Optical Telescope (NOT), and checked the 
images visually to locate new giant arcs. Thus the SGAS probes the comoving volume of 
~ 6 Gpc 3 , representing the largest giant arc survey conducted to date. Although the survey 
is still on going, it has alr e ady u ncovered more than 30 new lensing clusters that exhibit giant 
arcs. See iHennawi et al.l (120081 ) for more details on the survey method and initial results of 
the SGAS. 

In this paper, we study the following four clusters; Abell 1703 (hereafter A1703, z = 
0.281), SDSS J1446+3032 (SDSS1446, z = 0.464), SDSS J1531+3414 (SDSS1531, z = 
0.335), and SDSS J2111— 0115 ( SDSS2111, z = 0.637) . Giant arcs of all the fo u r clusters were 
newly disco vered by our su r vey (Hennawi et al.l 120081 ). iLimousin et al.l (120081 ). ISaha &: Read 
(120091 ). and iRichard et al.l (120091) conducted strong lens ing analysis of A1703 based on the 
Hubble Space Telescope images. iBroadhurst et al.1 (120081 ) presented results of mass modeling 
of A1703 from both strong and weak lensing. Besides these, no lensing analysis has been 
published for these clusters. 



2.1. Imaging Follow-up 



We conducted imaging observations of these clusters with the Suprime-Cam ( IMiyazaki et al. 



20021 ) on the Subaru 8.2-meter telescope on 2007 June 15. It covers a field of view ~ 34' x 27' 



with the pixel scale of 0'/202. For each cluster, we obtained deep g-, r-, and z-band images. 
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We assign longer exposure time to r-band images because we conduct weak lensing analysis 
using r-band images. The g- and z-band images are taken to study colors of strong lens 
candidates and to select galaxies for shear measurements. We summarize the exposure time, 
the seeing size, and the li miting magnitude of each image in Table [TJ The imaging data are 
reduced using SDFRED (jYagi et al.ll2002l ; lOuchi et al.ll2004l ). The photometric calibration 
is performed using the SDSS data by comparing magnitudes of stars between SDSS and 
Subaru image s. Astrometri c calibration is performed using reference objects in the USNO- 
A2.0 catalog (IMonetl Il998t ). The resulting astrom etric accuracy is ~ O'M . We construct 
object catalogs for these images using SExtractor (IBertin fc Arnoutslll996l ). In this paper 
we basically use total magnitudes (MAG_AUTO). However, for measurements of galaxies 
colors ( §4.2p we adopt aperture magnitudes (MAG_APER) with an aperture diameter of 
2", because aperture magnitudes tend to provid e better measuremen ts of colors. All the 
magnitudes are corrected for Galactic extinction (jSchlegel et al.lll998l ). 



2.2. Spectroscopic Follow-up 



Spectroscopic observations were conducted for three of the four clusters studied in this 
paper. SDSS1446, SDSS1531, and SDSS2111 were observed with the Frederick C. Gillett 
Telescope (Gemini North) between the months of 2008 February and 2008 July. The primary 
goal of the spectroscopic observations was to obtain redshifts of arcs to facilitate strong 
lensing modeling. We briefly summarize the spectroscopy here, and refer the reader to 
Hennawi et al. (2009, in preparation) for additional details. 

All spectro scopic observation s were carried out using the Gemini Multi-Object Spectro- 
graph (GMOS; iHook et al.l 120041) using multi-object sl i tmasks in microslit n od- and- shuffle 
(N&S) mode (e.g jGlazebrook fc Bland-Hawthornll200ll ; lAbraham et al.ll2004j ). Custom slit- 
masks were designed and arcs were targeted based on their colors in the deep Subaru imaging. 
After targeting all of the arc candidates, any remaining slits were placed on cluster members, 
easily identified by their red sequence colors. 

Spectra were taken with the R150_G5306 grating in first order which gives a dispersion 
of 3.5 A per pixel, with six pixels per resolution element resulting in a spectral FWHM~ 
940kms _1 . Although the R150 grating offers broad spectral range from the atmospheric 
cutoff to A > 1/mi, the drop in sensitivity at the blue and red extremes, due both to 
the GMOS CCD and the R150 grating efficiency, results in effective spectral coverage of 
4000 - 9500A. Our exposure times were 2400 sec resulting in 1200 sec effective integration 
for each of the two submasks. Three exposures were taken for each target. Thus if an arc 
was targeted on both submasks (typical for the most promiminent arcs) the total integration 
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time was 7200 sec. 



3. Strong Lensing Analysis 



3.1. Data 



We obtain strong lensing constraints of the clusters by fitting the positions of arcs. 
Figure [TJ presents the Subaru S uprime-cam images of the central regions of the clusters. As 
shown in lHennawi et al.l (120081 ). these clusters exhibit clear giant arcs. We identify multiple 
images based on their colors, redshifts, and also iteratively by building preliminary mass 
models of the clust ers. For A1703, we ad opt the identification and spectroscopic redshifts 
of arcs reported by iLimousin et al.l (120081 ) and iRichard et al.l (120091 ) . The redshifts of giant 
arcs of SDSS1531 were successfully measured with the Gemini telescope. Although we could 
not measure the redshift of the lensed blue images in SDSS1446 from our Gemini spectrum, 
the absence of [Oil] 3727A emission line (which should be detected for blue, actively star- 
forming galaxies like this) and the brightness of the arc in the g-band image suggests that 
the redshift should be in the range 1.6 < z < 3.5 (Hennawi et al. 2009, in preparation). 
We adopt this redshift range in our mass modeling. Since the spectrum of the giant arc of 
SDSS2111 was inconclusive, we conservatively assume its redshift to be z < 3.5 based on 
the g — r color. Although some additional possible multiply imaged systems are identified in 
some of these clusters, we conservatively restrict our analysis to those image systems where 
lensing hypothesis is completely unambiguous. We summarize the locations and redshifts of 
the multiply imaged systems for strong lens modeling in Table |2J 



3.2. Mass Modeling 

We model the dark matter distribution in each clusters by the NFW profile: 

{r/r s ){l + r/r s ) z 

where p s is a characteristic density and r s is a scale radius. Throughout the paper, we adopt 
the virial overdensity A v ; r (;2:) to compute the virial mass M v [ r and virial radius r v i r : 



4 



vir 



M vir = jnrl iv A viT (z)p(z) = / 4vrr 2 p(r)cir, (2) 



o 



where p(z) is the mean matter density of the universe at redshift z. The nonlinear overdensity 
A v i r (2;) is calculated at each z in a standard way using the spherical collapse model. The 
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concentration parameter c v ; r of the model is defined by 



r 



vir 



(3) 



Thus clusters with larger values of c v ; r are more centrally concentrated. iV-body simulations 
based on the ACDM model predict that concentration parameters depend on halo masses 



and redshifts; for massive clusters, th ey are typica 



redshift as (1 + z)~ a with a ~ 0.7 - 1 (IBuI 



20031 : iHennawi et aUBOOtt iNeto et al.ll2007l : iDuffv et al. 



ock et al. 



y ~ 5 at z = and evo 



2001 



Wechsler et al 



2008 



2002 



Maccio' et al. 



ve with the 



Zhao et al. 



2008|). 



To compute lensing properties, we first project the spherical NFW profile (eq. pQ) 
onto a two-dimensional lens plane. Next we include an ellipticity e by replacing the radius 
r = a/ x 2 + y 2 in the projected mass density to — e)x 2 + y 2 / (1 — e). Then we allow 
an arbitrary position angle 9 e by rotating the projected mass density. T he lensing deflec- 
tion angle for our mass profile is computed using the method described in ISchramml ( 1990 ). 
Since the ellipticity is introduced in the surface mass density rather than the lens potential, 
our model does not suffer from unphysical mass distributions (such as dumbell-like isoden- 
sity contours a nd negative mass den sities), which may be seen in the case of elliptical lens 
potentials (e.g.. lGolse fc Kneibll2002j ). 



We also include member galaxies in our mas s modeling. The total mass distribution 
of each member galaxy is modeled by the pseudo- lJaffd (119831 ) model which has the three- 



dimensional radial profile of p oc r~ 2 (r 2 + r^J" 1 , i.e., the isothermal model with the trunca- 
tion at r cut . We identify red member galaxies in the color-magnitude (g — i versus r) diagram 
and include ~ 40 brightest member galaxies for the mass model of each cluster. We fix the 
locations, ellipticities, and position angles of the members to those measured in the Subaru 
Suprime-cam data. Note that including the central galaxies is important for accurate esti- 
mates of dark matter core masses, as their effects on strong lensing c an be significant (e.g., 
Meneghetti et al.ll2003l ; IHennawi et al.l 120071 ; IWambsganss et al.l 120081 ) . On the other hand, 
the effect of other member galaxies is rather minor, but they can be important in accurately 
reproducing positions of some lensed images. 

To reduce the number of parameters, we adopt the exact sca ling relations for the v elocit y 
dispersions a and cutoff radii r cut of member galaxies following iNataraj an &: Kneibl (119971 ). 
who constrained the mass distribution around cluster elliptical galaxies using detailed strong 
and weak lensing models. Specifically, we assume that they scale as a oc L 1//4 (Faber- 
Jackson relation) and r cut oc L 1 ! 2 (corresponding to a constant mass-to-light ratio), where 
L is a total r-band luminosity of the galaxy. We apply this scaling relation to all member 
galaxies including the BCG of the cluster. To avoid unrealistically massive member galaxies, 
we add a Gaussian prior to the normalization of the scaling for the velocity dispersion, 
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which is estimated using the correlation with the magnitude (IBernardi et al.l 120031 ). but 
we adopt a conservative 20% error so that a wide range of mass models are allowed. We 
also estimate the normalization of the scaling for the cutoff radius by that inferred from 
observations, r cut = 18(cr/200kms -1 ) 2 k pc (where we fix a in this relation to the mean 
normalization of the Gaussian prior on a) (INatarajan et al.ll2002l . |2009| ). but again allow r cut 
to be Gaussian-distributed with a 50% error which crudely models the uncertainty of this 
relation in observations. 

Finding multiple images for a given source position requires intensive search in the 
image plane. This is time-consuming, even if we adopt adaptive-mesh refinement scheme to 
improve grid resolutions only near the critical curves, preventing us from exploring a large 
parameter space in a reasona ble time scale. I nstead, in this paper we evaluate \ 2 i n the 
source plane to speed up (e.g.. lKochanek!ll99ll ): 



y 2 
Asrc 



\(JLi(Ui - U s 



a 



+ prior, 



(4) 



pos 



where /ij is the magnification tensor for observed i-th image, u; are the source position 
computed from the position o f the i-th ima ge, and u src is the best-fit source position which 
can be found analytically (see lKeetonll200ll ). Since the magnification tensor is nothing but a 
mapping from the source plane to the image plane, |/ii(uj — u src )| approximates the distance 
between observed and predicted image positions in the image plane; indeed, we confirmed 
that this source plane \ 2 is reasonably accurate at around our best-fit models, compared 
with the standard image plane x 2 . Note that this approach automatically ignores any extra 
images. An important parameter here is the positional uncertainty in the image plane, cr pos , 
which is directly related with errors on best-fit model parameters. While the measurement 
error was sufficiently small, comparable to the pixel scale of the image (~ 0'.'2), it has been 
known that m ultiple images can be fitt ed only with ~ 1" accuracy for some massive lensing 



clusters (e.g., iBroadhurst et al.l l2005bl ) . presumably because of complex nature of cluster 



mass distributions. Thus, in this paper we assume the positional error of 1'.'2 (image plane) 
for all multiple images. 

We explore the mult i- dimensional x 2 surface using a Markov chain Monte Carlo (MCMC) 
approach. We adopt a standard Metropolis-Hastings sampling with the multivariate- Gaussian 
as a proposal distribution. During the sampling, the range of the virial mass is restricted to 
1.4 x 1O 14 M < M vir < 7 x 1O 15 M , and that of the concentration parameter to c vir < 40. 
We then derive constraints on parameters by projecting the likelihood distributions to the 
parameter space. Mass modeling is performed using the software glafic (M. Oguri, in prepa- 
ration). 
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3.3. Result of Strong Lens Modeling 



Results of our strong lens modeling are summarized in Table [31 and the critical curves 
of the best-fit models are displayed in Figure HJ First, we find that our best-fit models 
successfully reproduce the multiple images with reduced chi-squares of x 2 /dof < 1, which 
suggests that our assumption of the positional error of l'/2 was reasonable. Second, strong 
lensing constrains the center of the dark halo component quite well with a typical error of an 
arcsecond, and the best-fit dark halo center is close to the position of the BCG. An exception 
is SDSS2111 for which the offset of the center of the dark halo from the BCG is large, > 10" 
(see below). Third, it is found that the virial mass M vir and concentration parameter c v ; r 
are not constrained very well. In fact this is expected, because multiple images of a source 
redshift mainly constrain the mass enclosed by the arcs, which is sensitive both to M v j r and 
c v ; r . Fourth, we find that relatively large ellipticities of e ~ 0.4 are required to fit the images. 
An exception is SDSS1531 for which the mass distribution appears to be rather circular, 
e ~ 0.1, although an ellipticity as high as ~ 0.5 is a llowed. These large ellipt icities are in 



fact broadly consistent with ACDM predictions (e.g.. lOguri fc Blandfordll2009l ) 



Our best-fit mass model of A1703 from stro ng lensing (see also Appendix IA1 a 



to be in reasonable agreement with earlier work (jLimousin et al.l 120081 ; iSaha fc Read! 12009 



ap 
3 



pears 



Richard et al.l 120091 ) . For example, best-fit ellipticity and positio n angle of dar k halo com- 
ponent, e = 0.36 and 9 e = —24.1, agree well with the result of iRichard et al. J2009h. the 



ellipt icity in the lens potential of e pot = 0.23 and 9 e = —26.0. (see, e.g., iGolse &: Kneib 



2002 



for the relation between e pot and e). Moreover, we find that the structure of the 
critical curve of A1703 is consistent with those found in the earlier work. It is worth noting 
that the central ri ng (A1-A4) represents a ra re exa mple of lensing by a hyperbo l ic um bilic 



catastrophe (e.g.. iBlandford fc Narayanlll986l ). See lOrban de Xivry fc Marshalll (120091 ) for 



more comprehensive discussions about the central ring as a hyperbolic umbilic catastrophe. 

One of the most important quantities that are well constrained by strong lensing is 
the Einstein radius, #Ein- We obtain constraints on 6*Ei n at the redshifts of the arcs from 
our strong lens modeling as follows. For all the MCMC samples we estimate the Einstein 
radii of the dark halo (NFW) components by forcing e = and calculating the radii of the 
outer critical curves predicted by the dark halo components aloneJl] Thus the contribution 
of cluster member galaxies including the central galaxy is removed in our estimate of 



*We compared radial profiles of dark halo components with e = with azimuthally averaged radial 
profiles of elliptical dark halo components, and found that they agree with each other at <3% level, which is 
negligibly small compared with statistical errors. Thus the procedure adopted in this paper does not involve 
any significant systematics. 
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suggesting that the value could in fact be smaller than observed distances of lensed images 
from the center. We emphasize that this quantity serves as a useful strong lens constraint 
when combined with weak lensing results, particularly because we are interested in the radial 
profile of the dark matter component, rather than the total matter including stars which 
dominate at the very center of clusters (see $5])E| The derived constraints on 6*Ei n are shown 
in TableHl Since multiple strong lens systems with different redshifts are available for A1703, 
for this particular cluster we derive Einstein radii at two different source redshifts, z = 0.8889 
(the redshift of arc Al— A5) and 2.627 (the redshift of arc Bl— B3 and Cl— C3). We regard 
these two Einstein radii as statistically independent, which is reasonable given the large 
difference of 9 Ein between these two redshifts. 

Here we briefly discuss the mass model of SDSS2111, which appears to exhibit the 
large offset between the lens potential center and the location of the BCG (see Figure (TJ 
the best-fit lens potential center corresponds to the center of the inner critical curve). We 
find that no acceptable lens model has the offset smaller than ~ 5". The offset is required 
because the center of curvature of the long bright arc is clearly displaced from the BCG. 
However, as seen in Table El the current data do not constrain the center of the dark halo 
(NFW) component very well, which is the main reason that #Em involve relatively large 
errors. Thus it is important to improve the strong lens constraints by adding more multiply 
imaged systems. 



4. Weak Lensing Analysis 
4.1. Distortion Measurements 



We derive weak lensing shear signals following the algorithm outlined in iKaiser et al. 



(119951 . hereafter KSB). First, for each cluster we reprocess Subaru Suprime-cam r-band 



images through IMCAll to compute shapes of each objects. In order to compute shapes 
accurately, we iteratively refine the centroid of each object. During the process objects with 
offsets in each iteration larger than 3 pixel are removed. The ellipticities e of objects are 
measured from the weighted quadruple moments of the surface brightnesses. We select stars 
for Point Spread Function (PSF) correction in a standard way by identifying the stellar 
branch in the magnitude versus half light radius plane. We apply an additional cut, 



2 Strictly speaking, the Einstein radius derived in this way actually constrains the mass distribution of 
dark matter plus hot gas, rather than the dark matter distribution alone. 

3 ht tp : / / www . ifa. hawaii . edu / ~ kaiser /imcat / 
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the signal-to-noise ratio v > 15 and the ellipticity |e| < 0.1, to construct a final sample of 
template stars for the PSF correction. We correct ellipticities of objects using the sample of 
stars as reference 



e - P sm q* 



(5) 



where P sm is the smear polarizability tensor and q* = (P s * m ) _1 e* is the PSF anisotropy kernel 
estimated from the template s tars (hereafter the supers cript * denotes quantities m e asure d 
for template stars). Following lOkabe fc Umetsul (120081 ) and lUmetsu fc Broadhurstl (120081 ). 
we obtain a smooth map of e* by first dividing the whole image into 5x4 chunks and fit e* 
in each chunk using second-order bi-polynomials with iterative 3cx-clipping rejections. We 
show ellipticities of the template stars before and after the PSF correction in Figure [2] and 
Table [5) The Figure demonstrates that our procedure successfully reduces PSF anisotropies 
over the entire field. Residual stellar ellipticities after the PSF correcti on are comparable to 



those in other weak l ensing studies using Subaru Suprime-cam images (IHamana et al.l 12003 



Mivazaki et~al1l2007l : bkabe fc Umetsul I2OO8I : lUmetsu k Broadhurstl [io08h . 



Next we correct the isotropic smearing effect to estimate weak lensing shear. The pre- 
seeing shear polarizability P g relates observed galaxy ellipticities with their true values, and 
is calculated as 

(6) 



Pg — Psh Psm(P sm ) P s h; 



where P s h is the shear polarizability tensor. Since the measurement of P„ for each galaxy is 



extremely noisy, we appr oximate P g as P g times an identity matrix with (jErben et al.l 12001 
Hetterscheidt et al.ll20Q7h : 



P! 



2 i tr ( p shj 



tr(P< 



^r(P sl 



(7) 



Here we use a smooth map of tr(P s * h )/tr(P s * m ) constructed by fitting them as a function of 
the position, with a fifth-order polynomial over the entire template stars. We then derive 
the reduced shear g = 7/(1 — k) for each object as 



cal " 



We use raw, u nsmoothed P! bec ause the smoothing does not necessarily improve the shear 
estimate (e.g.. lErben et al.ll200ll ). We include the calibration correction factor / cal = 1/0.88 
since the KSB algorithm with t he scalar correction scheme for P g has been known to under - 
estimate shears by ~ 10 — 15% (jErben et al.ll200ll ; iHeymans et al.ll2006l : iMassey et al.ll20071 ). 
We confirmed that our shear estimate procedure, together with the source galaxy selection 
described in §4.21 recovers lensing shears in the simulated Su baru Suprime-cam images of 
the Shear Testing Programme 2 (STEP2; IMassey et al.l 120071 ) quite well, with typically a 
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few percent of the shear-calibration bias (parameter m) and an order of magnitude smaller 
residual shear offset (parameter c). 



4.2. Source Galaxy Population 

The background galaxy selection is crucial for weak lensing analysis. We choose following 
galaxies for weak lensing analysis; (1) no close companions at distances within 10 pixels, (2) 
the half light radius larger than 1.2r£, where r* h is the median of the half light radii of 
template stars used for the PSF correction, (3) the half-light radius smaller than 10 pixels, 
(4) the signal-to-noise ratio v > 7, (5) the pre-seeing shear polarizability PJ > 0.05, and (6) 
the shear value \ j\ < 2.0. In addition, we apply magnitude and color cuts, as we will discuss 
below. 

An appropriate c olor selection is import ant to minimize the dilution effect by cluster 



member galaxies (e.g.. iMedezinski et al.ll2007l ). One way to avoid contaminations from clus- 



ter members is to use galaxies redder than the red sequence of the cluster or bluest galaxies. 
In this paper we adopt blue galaxies, because our clusters are located at moderately high 
redshifts (z ~ 0.3 — 0.7) and hence very few red background galaxies are available. Specif- 
ically, we use galaxies with g — i < 1.0 (measured in aperture magnitudes) for our weak 
lensing analysis. This is much bluer than red-sequence galaxies (g — i~ 1.7 — 2.8 for r ~ 24), 
suggesting that the dilution effect is not significant for our source galaxy samples. 

We now estimate the depths of our galaxy sample for weak lensing analysis. A key 
quantity here is the mean lensing depth averaged over the population of source galaxies: 

dis \ _ f dz d Pw\ dis / g N 
dosl J dz d os } 

with di s and d os being angular diameter distances from the lens to the source and from the 
observer to the source, respectively. The PDF dp w \/dz denotes the redshift distribution of 
the source galaxy sample, which is particularly important in our study because of relatively 
high redshifts of the clusters. 

To estimate dp w \/dz, we adopt a photometric redshift sample of galaxies in the C anada- 



France Hawaii Telescope Legacy Survey (CFHTLS) presented by lllbert et al.l (120061 ). The 
idea is that we apply our color and magnitude cuts to the CFHTLS galaxy sample, which 
is well calibrated by the large amount of photometric and spectroscopic data, to infer the 
redshift distribution in our Subaru galaxy sample. First, we compare the number counts of 
galaxies with g—i < 1.0 between the CFHTLS and our source galaxy samples. Figure [3] shows 
the relative weight as a function of total r-band magnitudes, w(r) oc ^subam(0/ n CFHTLs( r )) 
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i.e., the number ratio of the source galaxies and CFHTLS galaxies. Sharp declines at r ~ 25 
imply that the effect of our background galaxy selections is significant beyond this magnitude. 
We therefore use galaxies with total r-band magnitudes 22 < r < r max for our weak lensing 
analysis, where r max is chosen so that the relative weight becomes approximately 0.3 (see 
Table [6] for exact values). Next we estimate the probability distribution functions (PDFs) of 
redshifts for our background samples dp w \/dz as 



dz 



Z>( r ) 



(10) 



where the summations run over all CFHTLS galaxies with g — i< 1.0 and 22 < r < r ma , x , and 



dp/d z(r) is the PDFs of photometric redshifts for individual CFHTLS galaxies (jllbert et al. 
20061 ). With this PDF, we estimate the lensing depths (di s /d os ) from equation (Q. 



Table M summarizes the mean depths of our source galaxy sample. We also compute 
the effective source redshift z<i that is defined by the redshift at which the distance ratio 
di s /d os becomes equal to the mean distance ratio {di s /d os ); in most situations of weak lensing 
analyses we can simply assume that all the source galaxies lie at z — Zd- However, since 
the reduced shear g = 7/(1 — n) is not exactly proportional to di.J d nx , adopting a common 
redshift for all backgr ound galaxies in fact underpredicts values of g (ISeitz fc Schneiderlll997l ; 
Hoekstra et al.ll2000l ). Thus, in computing expected g for a given mass model we fully take 
account of the redshift distribution of background galaxies dp w \/dz instead of adopting the 
effective source redshift z&. 

9000 — 10000 (surface number 
2 ), except SDSS1531 for which the number of source galaxies is 
7 arcmin -2 ), because of the worse seeing size of the image. The 



The resulting numbers of background galaxies are 
density of n g ~ 10 arcmin 
much smaller, ~ 6000 

reason for the relatively small number density compar ed with the expected number density 



\P>g ~ 



of n g ~ 30 — 40 arcmin 2 for Subaru weak lensing (e.g.. iMiyazaki et al.ll2007l ) is our color cut 



to eliminate the dilution effect. In fact, the number density is com parable to other Subar u 



cluster weak lensing studies which adopted similar color cuts (e.g.. lOkabe fc Umetsull2008l ). 



4.3. Two-Dimensional Mass Distributions 

Before studying radial profiles, we check two-dimensional mass (conv ergence) distribu- 



tions r econstructed from weak lensing shears using an inversion algorithm of lKaiser fc Squires 



( 119931 ). Here we assume weak lensing limit (k< 1, i.e., g ~ 7), which suggests that the 
reconstructed k will be slightly overestimated in the cores of clusters. We use these two- 
dimensional mass maps only to study the morphology of the cluster; we constrain cluster 
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mass distributions using radial profiles of tangential shears (see §4.4p . We also derive lumi- 
nosity density maps of red-sequence galaxies for comparison. 

We show Gaussian-smoothed (a = 60") maps of all four clusters in Figure HI In all 
the clusters, weak lensing signals are clearly detected. In addition, we find that the recon- 
structed mass distributions well trace the spatial distributions of red-sequence galaxies, with 
roughly similar mass and light centroids. The clusters exhibit deviations from circular mass 
distributions, in particular for A1703. The umidordal structure with aligned position angles 
of weak lensing, the light distribution, and strong lensing (see the critical curve shown in 
Figured]) implies that the cluster A1703 is quite relaxed. For SDSS1531, we find a probable 
secondary peak at the ~ 4' northwest of the cluster center, implying that the cluster may 
not yet be relaxed. The position angles of mass distributions from strong and weak lensing 
appear to agree. Several subpeaks and non-regular morphology of SDSS2111 suggest that 
this cluster is not relaxed either. However, SDSS2111 again shows agreement between weak 
and strong lensing in that the isodensity contours appear to be elongated in the north-south 
direction. In addition, the centroid of the mass distribution is slightly offset North of the 
BCG, which is consistent with the result of strong lens modeling. It is important to obtain 
X-ray data for these clusters to explore their dynamical states further. 



For each cluster we derive an azimuthally averaged one-dimensional shear profile to 
obtain constraints on the cluster mass profiles. The center of each cluster is fixed to the 
center of the dark halo component determined from strong lens modeling (see Table [3]). For 
a given center, we can calculate the tangential shear g + and the 45° rotated component g x 
from the reduced shear g = (gx, g 2 ) as follows: 



where is the polar angle. We obtain g + and g x for each radial bin by simply averaging 
these shears of individual source galaxies without weighting. We estimate the mean and 
error using a jackknife resampling technique. The shear profiles are computed in the range 
of radii 9 = [0'.8, 20'] with a bin size of Alog# = 0.2. Thus we have seven radial bins in 
total. 

In Figure [5] we plot radial profiles of g + and g x of all the four clusters. We detect 
weak lensing signals significantly up to ~ 10' from the cluster center. On the other hand, g x 
is consistent with a null signal, as expected for shears produced by gravitational lensing of 



4.4. Tangential Shear Profiles 



9+ = -pi cos 20 - #2 sin 20, 
9x = -9i sin 20 + g 2 cos 20, 



(11) 
(12) 
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clusters. We fit the tangential shear profile by the NFW profile (eq. [I]) with the virial mass 
M vir and concentration parameters c v ; r as free parameters. For this, we adopt the following 

X 2 : 

2 _ \- - 9+(0i\ M vir, C vir )] 2 , , 

XWL — / , 2 ' y 16 ) 

*-r* erf 

where and exj denote observed tangential shear and its error for z-th radial bin. As 
shown in Figure [5j the NFW model provides reasonable fits to the data. The best-fit model 
parameters are summarized in Table [7J 



5. Combining Strong and Weak Lensing 

5.1. Constraints on Parameters 

In this section, we combine results from strong lensing ($3]) and weak lensing (£j4j) to 
constrain the radial mass profiles over a wide range of radii. We do so simply by summing 
up the chi-squares: 

X 2 = Xk + Xwl, ( 14 ) 
where the weak lensing constraint Xwl * s fr° m equation ( fT3l) . Although we can in principle 
use strong lens constraints obtained from detailed fitting of image positions in here 
we choose a conservative approach to adopt only the Einstein radius as our strong lens 
constraint. Specifically, Xsl * s computed as 

2 _ [$Ein — #Ein(^vir, 

XSL = 2_^ ' 

where ^Ein and o"Ein are the Einstein radius and its errors listed in Table HI We include the 
asymmetry of errors on ^ Ein by using different values of between 6 > #Ein and 6 < ^Ein- 
For SDSS1446 and SDSS2111, only the allowed range of the redshift of the arcs (and hence 
the redshift of the Einstein radius, ^Em) was given; therefore, for these clusters we assume a 
flat prior 1.6 < z^m < 3.5 (SDSS1446) or z Ein < 3.5 (SDSS2111) and marginalize x|l over 
2Ein to derive strong lens constraints. We note that the contribution of stars (member galaxies 
and the BCGs), which can be significant at strong lensing regime, is excluded in estimating 
^Ein from strong lens modeling. Thus the derived radial profile should be considered as that 
for the sum of dark matter and intracluster gas rather than the total matter profile. 

We show the constraints in the M vir -c v i r plane in Figure [61 As expected, adding strong 
lensing significantly improve constraints on these parameters, especially the concentration 
parameter c v ; r , although this is not the case for SDSS2111 because of the unknown arc 



-16- 



redshift and the relatively large error on ^Ein- Constraints from strong and weak lensing are 
consistent with each other, except A1703 for which the best-fit models slightly differ. We 
list the best-fit model parameters in Table UJ It is found that the radial profiles are fitted 
well by the NFW profile with x 2 /dof ~ 1. 

To see how strong and weak lensing probe different radii, we plot radial profiles of con- 
vergence k from lensing observations and best-fit models in Figure UJ We derive convergence 
at each radial bin from weak lensing shears by using the following relation: 

g+{0) [1 - k(0)] = ^ jT 9>K(e')d8' - k{6). (16) 

We need a boundary condition in order to solve this equation. We fix k in the outermost bin 
to the value computed from the NFW profile w i th be st-fit parameter values from weak lensing 



data alone. See, e.g., lUmetsu fc Broadhurstl (120081 ) for practical procedures to compute k 



from discrete g + data. As clearly shown in the Figure, strong lensing probes radii roughly 
an order of magnitude smaller than those constrained by weak lensing. Hence by combining 
strong and weak lensing we can constrain the density profile of the cluster over 2 dex in radius. 
We find that the model that fits weak lensing data of A1703 underpredicts convergence in 
the core of the cluster; the Einstein radius of A1703 implies more centrally concentrated 
profile of A1703 than expected from weak lensing data alone. For the other three clusters, 
convergence profiles from weak lensing alone and strong and weak lensing agree quite well. 
However, strong lensing data do narrow down the allowed range of radial profiles, as is clear 
from Figure El 

We now compare our r esult of A1703 with the combined strong and weak lensing analysis 



by iBroadhurst et al.l (120081 ) . We find that the best-fit virial masses are consistent with each 
other, but the concentration parameter of our best-fit model, c v i r = 6.51q; 7 , differs from their 
best-fit value, c vir = 9.9l^g at ~ 2a level. We ascribe the difference to the larger Einstein 
radius they assumed, ^in = 33" for z s = 2.8. Indeed, our best-fit model predicts 0Em = 22" 
for z q = 2.8, which appe ars more consistent with the locations of lensed arcs in A1703 (see 



also iRichard et al.l l2009h . 



5.2. Distribution of the Concentration Parameter 

One of our main interests lies in the possible excess of the concentration parameter 
c v ; r found among lensing clusters. Here we compare our results with theoretical predictions 
based on the ACDM model. 

For the theoretical prediction, we adopt results of iV-body simulations in the Wilkinson 
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Microwave Anisotropy Probe (WMAP) 5-year cosmology (jDuffy et al.ll2008l ). They derived 
the mean concentration of dark halos as a function of halo mass and redshift as 



7.85 



sim 



1 + z 



,0.71 



2.78 x 10 12 M m 



-0.081 



(17) 



The apparent concentration parameters for halos with given mass and redshift are signifi- 
cantly scattered, with approximately the log-normal distribution with o"i ogc ~ 0.14. However, 
the diversity of the halo population found in A-body simulations suggest that clusters with 
giant arcs should represent a quite biased population. This lensing bias arises from the fact 
that the lensing cross section is a sensitive function of the halo concentration. Therefore 
we expect projected mass distributions of lensing clusters to be more centrally concentrated 
than ordinary clusters, implying that the clusters are intrinsically more concentrated and/or 
their major axes are preferentially aligned with the line-of-sight, giv ing rise to app a rently 
more concentrated distributions of projected surface mass densities. iHennawi et all (120071 ) 
estimated the lensing cl usters have ~ 40% higher values of con c entrat ions compared with 
normal clusters (see also lFedeli et all 120071 ). lOguri fc Blandfordl (120091 ) focused on clusters 
with larger Einstein radii and argued that such clusters have ~ 40 — 60% higher concen- 
trations. Based on these discussions, in this paper we consider 50% enhancement of c vir 
due to the lensing bias. Note that the lensing bias include s the p rojection effect from the 
halo triaxiality. From the calculation of lOguri &: Blandfordl ( 120091 ). it is also found that the 
apparent concentration of the lensing clusters has slightly smaller scatter with <J\ ogc ~ 0.12. 

Figure [8] shows the distribution of c v ; r of the four clusters we study in this paper as well 
as several lensing clusters studied before. We consider only clusters whose radial profiles are 
well constrained from combined strong and weak lensing analyses. We find that all the four 
clusters have best-fit value of c vir higher than the theoretical expectations. We still see the 
excess even if we take the lensing bias into account, although the value of c v ; r for each cluster 
is marginally consistent with the theory if the error is taken into account. Put another way, 
the exces s is not so strong as claimed i n earlier work based on analysi s of different lensing 
clusters ( iBroadhurst et al.l l2005al . 120081 : IComerford Natarajanl 120071 ) . By combining all 
the 10 clusters, it is clear that our sample of lensing clusters has larger concentrations than 
expected from the M vir -c v i r relation predicted by the ACD M model. The concentrations are 
also higher than the relation determined observationally by IComerford fc Natarajanl (120071 ) 
for a sample of lensing and X-ray clusters. Since our sample cannot constrain the redshift 
and mass dependence of the mean concentration very well, we fix them to those in equation 
(fT7l) and fit the overall normalization to the data. We find 



c(fit) 



12.4 



,0.71 



10 15 M, 



-0.081 
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as our best-fit to the results of 10 lensing clusters. The data are inconsistent with the concen- 
tration parameter predicted in the ACDM by 7cr, even if we include the 50% enhancement 
to account for the lensing biasQ A simple average of theconcentration parameter (with the 
inverse of the measurement error as a weight) ignoring the mass and redshift dependence is 
(cvir) = 9.3 ± 2.6. 

We also check the excess as a function of the Einste in radius, which is a central qu antity 
to characterize the central structure of a cluster (e.g., iBroadhurst fc Barkanal 120081 ). The 
plot in Figure [S] shows that the excess of Cv ir depends slightly on the Einstein radius. We fit 
the data by a power-law and find 



r(fit) 



2.4 



sim 



#Ein 

35" 



0.41 



(19) 



where #Ein is the Einstein radius for the source redshift z s = 3. The hypothesis that this 
ratio does not depend on the Einstein radius is rejected at 99% confidence level. The weak 
dependence of the excess on the Einstein radius is reasonable in the sense that concentrated 
two- dimensional mass distributions are required for cl usters to produce large Einstein radii 



le.g. 



Broadhurst et all 12008 : lOeuri & Blandfordl 120091 ) . 



6. Summary and Discussion 

In this paper, we have studied the mass profiles of four clusters by combining strong 
lens constraints with weak lensing shear measurements. We have drawn our sample from 
giant arc clusters new ly discovered by the SGAS, the new arc survey using the SDSS data 



( jHennawi et al.ll2008l ). We take advantage of follow-up wide-field Subaru Suprime-cam im- 
ages to obtain weak lensing constraints out to nearly the virial radii of the clusters. The 
central densities of the clusters are determined well from the arcs with the redshifts spec- 
troscopically measured partly from our ongoing program with the Gemini telescope. The 
technique allows us to study the radial mass profile over 2 dex in radius, which is essential 
for reliable extractions of concentrations of the clusters. 

We have found that the radial profiles are fitted well by the NFW profile. We have deter- 
mined the virial mass M vir and concentration parameter c vir of the four clusters accurately 



4 The conclusion is unchanged even if we relax the lower limit of the arc redshift of SDSS1446 which 
originated from the absense of [Oil] emission line. This is because the lower limit of the arc redshift 
corresponds to the upper limit of the cluster core mass and hence to the upper limit of the concentration 
parameter. 
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from the combined strong and weak lensing analysis. We confirmed that strongly lensed 
background galaxies with measured redshifts indeed improve constraints on mass profiles, 
particularly c vir . We have found the values of c vir for our 4 clusters to be c vir ~ 8, except for 
SDSS2111 whose radial profile was not constrained very well because of insufficient strong 
lens information. The values are slightly higher than the ACDM predictions, even if we 
take account of the lensing bias that clusters with giant arcs are more centrally concentrated 
(intrinsically and/or apparently due to the projection effect) than normal clusters, although 
the excess is not as large as that claimed in earlier work. By combining all the 10 clusters 
with strong plus weak lensing analysis available, we confirm a 7a excess of the concentration 
parameter compared with the ACDM prediction. We find that the excess is dependent on 
the Einstein radius of the system such that clusters with larger Einstein radii show larger 
excess of the concentration parameters. 

There are several possible explanations for the excess of the concentration parameter. 



One such explanation is the effect of baryons (e.g., Kazant zidis et all 120041 ; iGnedin et al. 



2004 ; iPuchwein et al.ll2005l ; iLin et al.ll2006l ; iRozo et al.ll2008l ). In particular, the adiabatic 
contraction associated with the baryon cooling can enhance the core density of dark matter 
and hence can increase the concentration of clusters. Another possibility is that theoretical 
predictions are not so accurate. In particular, the probability distribution of the concentra- 
tion parameter for very massive halos (> 10 15 /i~ 1 M Q ) and its redshift evolution has not been 
studied very much in A^-body simulations, and thus our theoretical predictions inevitably 
rely on extrapolations from lower-mass halos. For example, it has been claimed that the red- 
shift evolution of th e concentrations of most massive halos may be different from that of less 
massive halos (e.g.. lZhao et al.ll2003l ). Thus it is important to perform many realizations of 
large box-size A^-body simulations to improve the accuracy of theoretical predictions at the 
high mass end. A more exotic interpretation is that clusters form earlier than expected from 
the AC DM model, as the con centration is known to correlate with the formation epoch of the 
cluster (IWechsler et al.ll2002l ). Such modification of the formation epoch can f or instance be 



realized by considering early dark energy or primordial non-Gaussianity ( e.g. 
20041 : ISadeh et alil2007l : ISadeh fc Rephaelill2008l : lOguri fc Blandfordll2009h . 



Mathis et al. 



This paper h as presented initial r esults of detailed studies for our unique sample of 
giant arc clusters (IHennawi et al.ll2008l ). To investigate the structure of lensing clusters in 
a more systematic and statistical manner, it is of great importance to extend this research 
by applying the technique we have developed in this paper to other lensing clusters. In 
addition, it is important to improve constraints on the mass models of individual clusters by 
adding more data. Measuring redshifts of more arcs will significantly refine our strong lens 
modeling. In addition, dynamical information from the velocity dispersion measurement, as 
well as X-ray and Sunyaev-Zel'dovich signals, provide an important cross check of our mass 



models (e.g.. lMahdavi et al.ll2007l ; lLemze et al.ll2008l . 120091 ). We are planning these follow-up 
observations for our sample of giant arc clusters. 
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A. Strong Lensing Modeling of A1703 by a Generalized NFW Profile 



The multiple strong lens systems in A1703 for a wide range of source redshifts constrain 
enclosing masses at different ra dii. Thus we expect t he dark matter density profile of this 



cluster can well be constrained (ILimousin et al 



20081 ) . By using the technique described in 



§21 we re-perform strong lens modeling of A1703 with the dark matter distributio n replaced 
from equation ([!]) to a so-called generalized NFW profile (e.g.. |Jing fc Sutoll2000l ): 



p(r) 



(r/r s ) a (l + r/r s ) 3 - 



(Al) 



where a = 1 represent the original NFW profile. In Figure [HI we show the probability 
distribution of a obtained from mass modeling. The result, a = 0.9^1, is consistent with 
NFW, suggesting the validity of our assumption of a = 1 in our strong lens modeling in 
§[3] An inner slope as steep as a = 1.5 is clearly rejected. The resul t is quite consistent 
with the recent strong lens modeling of A1703 by iRichard et al.l ( 120091 ). a = 0.92^^1, but 
wi th much larger error . It is not clear why our constraint on a is much weaker than that 
of IRichard et al.l (120091 ); one of the reasons may be that we are conservatively using only 
spectroscopically confirmed, robustly identified multiple images for our mass modeling. We 
note that an inn er slope slightly s hallower than NFW was suggested in some other lensing 
clusters as well rtSand et al.l lioosl ). The asymptoptic inner slope s hallower than —1 has 



also been implied by recent high-reso lution A-body simulations (e.g.. iFukushige et al.ll2004 
Graham et Zlbood : bao et all [20081 ) . 
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Table 1. Summary of Subaru Suprime-cam images 



Name 


Band 


Exposure 


Seeing 


TOl im 






cor* 
1 Del J 


1 dl UoCOJ 




A1703 


9 


400x3 


0.96 


26.39 




r 


300x7 


0.80 


26.31 




i 


240x5 


0.86 


25.58 


SDSS1446 


9 


400x3 


0.84 


26.22 




r 


300x7 


0.82 


26.24 




i 


240x5 


0.92 


25.47 


SDSS1531 


9 


400x3 


0.90 


26.28 




r 


300x5 


0.98 


26.19 




i 


240x5 


1.00 


25.38 


SDSS2111 


9 


360x4 


0.82 


26.13 




r 


300x8 


0.60 


26.05 




i 


240x7 


0.52 


25.48 



Note. - The seeing indicates the FWHM in 
the co-added mosaic image. The magnitude limit 
mii m refers to the 5a detection limit of point sources 
for 2" aperture diameter. In stacking SDSS1531 r- 
band and SDSS2111 r-band images, we removed 
a few frames with bad seeing values; these are ex- 
cluded from the numbers of frames indicated in this 
Table. 
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Table 2. Multiple images for strong lens modeling 



Name 


Image 


Ax 


Ay 


Ar 


Redshift 






[arcsec] 


[arcsec] 


[arcsec] 




A 1 7DQ 
Al I UO 


A 1 


1 1 7 
— 1 1 . ( 


P a 

O.O 


1 1 K 
lo.D 


n 888Q 

u.oooy 




A 
AZ 


P Q 
— O.O 


8 A 
o.4 


in 1 
1U.4 






A 1 

AO 


— 4.0 


9 

Z.O 


K 1 
O.O 






AzL 
ri4 


1 n ^ 

— 1U.0 


n 8 

U.o 


i n 7 

1U. ( 






A ^ 
AO 


99 1 
ZZ. 1 


1 A Q 

— 14. y 


OP. 7 
ZO. ( 






J31 


— oo.O 


1 1 

— 11.0 


K 

oo.o 


9 R97 
Z.0Z I 




R9 
JJZ 


9^ 7 

— ZD. I 


91 Q 

— z i .y 


11 8 
OO.o 






R^ 
J30 


91 1 
Zl.O 


98 K 
— Zo.O 


IK P 
00.0 






PI 
111 


Q9 Q 

— oZ.o 


i c q 
— 10. o 


IK 7 
OO. 1 


9 fi97 
Z.0Z I 






/I 

— OU.4 


1 a n 
— lo.U 


IK 1 
OO.O 






P3 
Lit) 


91 3 
Zl.O 


98 K 
— Zo.O 


00.0 






m 
jji 


i /i 1 

14. 1 


Q9 7 
OZ. 1 


oo.o 


1 QDS 

i.yuo 




Uz 


9c; n 
ZO.U 


9fi 7 
ZO. I 


IP 
OO.O 






JJo 


1 n /i 

— 1U.4 


17 8 

o ( .o 


9 

oy.z 








17 r 
— 1 ( .0 


Q9 9 

oz.z 


oo.o 


9 1PI~\ 
Z.OOU 




"F9 


P Q 

o.y 


ou.u 


OU.o 






Ji/o 


Q/l /I 
o4.4 


7 9 

( .z 


Q ^ 1 

oo. 1 






"F1 
r 1 


90 A 
ZU.4 


1 8 n 

— lo.U 


97 9 
z / .z 


9 IKK 
Z.000 




T?9 


ii n 
11. U 


1 Q 8 
lO.O 


17C 

1 ( .0 






r 


IP 1 
— OO. 1 


1 K K 
10.0 


oy.o 






F4 


-7.1 


-21.3 


22.5 




SDSS1446 


Al 


-11.9 


14.9 


19.1 


> 1.6 




A2 


13.2 


11.6 


17.6 






A3 


-4.9 


-10.7 


11.8 






A4 


18.9 


-5.3 


19.6 




SDSS1531 


Al 


-11.5 


9.7 


15.0 


1.096 




A2 


3.0 


-10.5 


10.9 






A3 


4.4 


-10.1 


11.0 






A4 


11.1 


-0.6 


11.1 






Bl 


-13.3 


5.7 


14.5 


1.095 
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Table 2— 


-(Inn t, inn pn 




Name 


Image 


Ax 


Ay 


Ar Redshift 






[arcsec] 


[arcsec] 


[arcsec] 




B2 


6.8 


3.8 


7.8 


SDSS2111 


Al 


10.0 


-5.1 


11.2 




A2 


6.1 


-8.0 


10.1 




A3 


-13.6 


-8.9 


16.3 




Bl 


8.6 


-6.2 


10.6 




B2 


8.1 


-7.1 


10.8 




B3 


-16.3 


-6.8 


17.7 




CI 


12.9 


-4.5 


13.7 




C2 


-0.1 


-11.1 


11.1 




C3 


-11.6 


-9.9 


15.3 



Note. - The same alphabet of images indicates that they 
are associated with the same source, inferred based on the colors 
and modeling. Ax and Ay denote the relative locations from 
the BCG of each cluster (see Figured]); the positive directions 
correspond to West and North, respectively. The distance from 
the BCG is shown as Ar. Redshifts are spectroscopic redshifts 
(see text for more details). The arc ide ntification and spectro- 
sco pic redshifts for A170 3 are taken from iLimousin et al.l (120081 ) 
and iRichard et all J2009h . 



Table 3. Result of strong lens mass modeling 



Name 


M vir 


x c 


Vc 


e 


e 


C v ir 


a 


r cut 


X 2 /dof 




[1O 15 M ] 


[arcsec] 


[arcsec] 




[deg] 




[kms- 1 ] 


[arcsec] 





A1703 


2 4l+ - 70 


-1.5 ±0.5 


0.71™ 


^fi+°- 04 

u -°°-0.03 


-24.1+g-? 


5.1 ±0.7 


337l 34 


14.9 ±3.6 


35.5/24 


SDSS1446 


u.oo_ 28 


1 1 +L4 

- - 1 - 1.6 




n 42 +0.44 


-36.3l|° 


11.7l? 8 2 3 


231 ± 38 


4 Q+l-9 
^•°-2.0 


0.0/0 


SDSS1531 


4 SQ+ 2 - 61 

^• oy -4.25 


n 4 +LS 

u -^-0.5 


-o.8±?:S 


i o +0 - 41 


-42.2i 3 7 ;? 


q n +13.0 


2911^ 


8.4+^ 


0.7/2 


SDSS2111 


U.dd_Q 27 


-3.9 ± 1.5 


i3.3±i:g 


4^+ - 27 


-168.2i| 7 


16.0± 24 ;? 


268^ 




1.5/6 



Note. — The parameters x c and y c denote the center of a dark halo (NFW) component relative to the location of 
the BCG; positive x c and y c indicate the offset to West and North directions, respectively. The position angle 9 e is 
East of North. The values of the velocity dispersion a and the cutoff radius r cut in this table are those of the BCG; 
other member galaxies have values obtained from the exact scaling relation a oc L 1 / 4 and r cut oc L 1 ' 2 . Note that wc 
have assumed the arc redshifts of SDSS1446 and SDSS2111 to be z = 1.6 when running the MCMC; some of the 
best-fit values (e.g., M vir ) is dependent on this assumption. 
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Table 4. Einstein radii from strong lens mass modeling 



Name #Ein [arcsec] z^i 



A1703 


12 1 +L6 


0.8889 




22.51}-? 


2.627 


SDSS1446 


16.51^8 


1.6-3.5 


SDSS1531 


q 3+2-9 


1.096 


SDSS2111 


18.61^ 


< 3.5 



Note. — The Einstein radius 0Em is 
that predicted by the best-fit dark halo 
component, i.e., contributions from 
member galaxies (including the BCGs) 
are excluded in deriving 6-^ [n . See text 
for more details. 



Table 5. Ellipticities of template stars for the PSF correction 



Name 


e* x 10 2 


e* x 10 2 


eUr x 10 2 


62, cor X 10 2 


-^star 


A1703 


+0.88 ± 1.11 


+0.56+1.22 


-0.02 + 0.49 


-0.01 + 0.43 


558 


SDSS1446 


-1.08 ±0.82 


-0.35 + 0.76 


-0.01+0.45 


-0.02 + 0.34 


598 


SDSS1531 


+0.79 ±0.94 


+0.14 + 0.85 


-0.01+0.45 


+0.00 + 0.39 


876 


SDSS2111 


-0.87 ± 1.82 


-1.60 + 1.22 


-0.01+0.70 


+0.03 ± 0.44 


1591 



Note. - Medians and standard deviations of ellipticities for template stars 
used to correct PSF anisotropies. e* and e* or are ellipticities before and after the 
PSF correction. The last column shows the number of the template stars. 
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Table 6. Source galaxy population for weak lensing analysis 



Name 




Magnitude 


(di s /d os ) 






A1703 


0.281 


22.0 < r < 25.7 


0.637 


0.899 


9155 


SDSS1446 


0.464 


22.0 < r < 25.3 


0.470 


0.999 


9985 


SDSS1531 


0.335 


22.0 < r < 24.9 


0.546 


0.832 


6032 


SDSS2111 


0.637 


22.0 < r < 26.0 


0.383 


1.181 


9721 



Note. - The range of r-band magnitudes, lensing depth 
(di s /d os ), effective source redshifts z^, and the numbers of source 
galaxies iV ga i are shown. See text for details. 



Table 7. Result from weak lensing and combined strong and weak lensing analysis 





Weak lensing 




Strong 


and weak lensing 


Name 


M vir [lO 15 M ] 


C-vir 


X 2 /dof 


M vir [lO 15 M ] c vir 


X 2 /dof 


A1703 
SDSS1446 
SDSS1531 
SDSS2111 


l - uo -0.50 

0.83^;1 
0.5933 


3.3 + \i 
9.ltl\ 4 
11.5±f 4 5 

14 1+25-9 
1 ^- i -9.5 


2.7/5 
6.3/5 
8.0/5 
7.5/5 


0.83l°;| 
66 +0 - 29 

U.UU_ 24 

92 +0 - 41 


6.5^;? 

8-311? 
7.9i?;° 
14.1^ 3 9 


7.9/7 
6.4/6 
8.1/6 
7.5/6 



Note. — Weak lensing constraints come from radial profiles of tangential shear g + (see 
Figure E]). The Einstein radii inferred from strong lens modeling are used as strong 
lens constraints (see Table HI). Note that we have restricted the range of the concentration 
parameter to c v ; r < 40. The virial radii of the best-fit models from strong and weak 
lensing analysis are 2.56, 1.91, 1.89, and 1.80 Mpc for A1703, SDSS1446, SDSS1531, and 
SDSS2111, respectively. 



-33- 




Fig. 1. — Subaru Suprime-cam images of 4 clusters we study in this paper. North is up and 
East is left. The size of each image is 107" x 80". Multiple images we adopt for strong lens 
modeling are named by an alphabet and a number (e.g., Al; see also Table [2]). Critical curves 
of our best-fit models at the redshifts of multiple images (z = 2.627 for A1703) are drawn by 
solid lines (see §3.31) . The brightest cluster galaxy for each cluster is marked by "BCG". The 
J2000.0 coordinates of the BCGs are 13:15:05.24+51:49:02.7 (A1703), 14:46:34.03+30:32:58.7 
(SDSS1446), 15:31:10.63+34:14:24.9 (SDSS1531), and 21:11:19.36-01:14:22.9 (SDSS2111). 
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10% ellipticity 

Fig. 2. — PSF anisotropy maps of the filed of SDSS2111 measured from template stars. 
The upper panel shows the raw ellipticities of the stars, whereas the lower panel is stellar 
ellipticities after the PSF correction is applied. The length and orientation of each line 
segment indicates the ellipticity and position angle of each template star. The size of each 
panel is ~ 35?5 x 2Q'.5, i.e., the entire field of the mosaic image is shown here. Table |5] lists 
ellipticities before and after the PSF correction for the fields of all the clusters. 



0.1 



A1703 
SDSS1446 
SDSS1531 
SDSS2111 




22 



23 



24 
r 



25 



26 



Fig. 3. — The relative weight w(r) of our background galaxy samples as a function of total r- 
band magnitudes, derived fro m the number rati o of source galaxies for weak lensing analysis 
and galaxies in the CFHTLS (jllbert et al.ll2006l ). The color cut g — % < 1.0 is applied to both 
galaxy samples. The weights are normalized by the maximum weight at r > 22. 



- 36 - 




[arcmin] [arcmin] 



Fig. 4. — Two-dimensional mass maps (solid) reconstructed from weak lensing shear mea- 
surements. North is up and East is left. Contours are drawn with spacing of Ak = 0.04, 
starting from k = 0.04. Note that the reconstructed k are slightly overestimated in the 
cores of clusters because we assumed the weak lensig limit, g ~ 7. The la noise level of 
each reconstructed mass map is a K ~ 0.03. Crosses indicate the locations of the BCGs. 
Contours of constant luminosity densities of red-sequence galaxies (dotted) are also shown 
for comparison. Both mass and luminosity maps are Gaussian-smoothed with a = 60". 
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10 
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Fig. 5. — Azimuthally-averaged tangential shear g + (eq. [IT] ) and the 45° rotated component 
g x (eq. [12]) as a function of distance from the cluster center. The NFW models fitted to 
the observed shear profiles are shown by solid lines (see also Table [7]). 
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M vir [MJ M vir [M ] 



Fig. 6. — Contours at la and 2a confidence levels in the M v ; r -Cvi r plane. Black and grey 
lines indicate constraints from weak or strong lensing, respectively. Combined strong and 
weak lensing constraints are plotted by shaded regions. The best-fit model parameters from 
combined strong and weak lensing are shown by crosses (see also Table [7]). 
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Fig. 7. — Radial profiles of convergence k for the source redshift z s = 3 constrained from 
strong and weak lensing. Radial profiles of best-fit models from combined strong and weak 
lensing are plotted by solid lines, whereas those from weak lensing only are shown by dashed 
lines (see also Table [7]). Filled squares are k profiles reconstructed from weak lensing tan- 
gential shear profiles shown in Figure [5j as boundary conditions, we assumed values of 
convergences in the outermost radial bin to be those computed from the best-fit NFW pro- 
files from weak lensing. We indicate the Einstein radii used as strong lens constraints by 
open squares, and the virial radii of the best-fit models from strong and weak lensing analysis 
by arrows. 
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20 30 40 50 
Ein( z s = 3 ) [arcsec] 



60 



Fig. 8. — Distribution of the concentration parameter c v j r constrained from combined 
strong and weak lensing analysis. Filled triangles show our r esults presented in this paper , 
wherea s filled squares are from literature, results for A1689 flUmetsu fc Broadhurstl 120081 ) . 
A2261 JUmetsu et al-l liooiih. A370, RXJ1347, CL0024 faroadhurst et al.ll2008h . and MS2137 
( IGavazzi et al. I l2003h . Upper: The distribution of 'unction of virial mass M vir . The 

data points are corrected for the redshift evolution of t he concentration p arameter predicted 
in the ACDM model, c vir oc (1 + z)~ a with a = 0.71 (fouffv et alihoOSl The grev shade d 
regions indicate la range of c v - ir derived from ACDM simulations by iDuffy et al.l (120081 ). 



Black shaded regions show predicted c oncentration parameters after approximately taking 
the lensing bias (~ 50%) into account (IHennawi et al.l 120071 ; lOguri fc Blandfordl 120091 ). The 
dashed line plots t he M v ] r -c v ] r relation determined observationally from a sample of lensing 
and X-ray clusters ( jComerford fc Natarajanl 120071 ). The dotted line shows the best-fit curve 
to the data. Lower: The ratios of c v ; r from analyses of lensing clusters to those expected 
in ACDM are plotted a function of the Einstein radius #Ei n for the source redshifts z s = 3. 
Again, the shaded regions show la range expected in ACDM, with (black) and without (grey) 
the lensing bias. The best-fit curve by a power-law is shown by the dotted line. 
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a 



Fig. 9. — Normalized probability distribution of the inner slope of the dark matter density 
profile, a (see eq. |Alj ). from strong lens modeling of A1703. Other model parameters are 
marginalized over. The vertical dotted line indicates the slope of the original NFW profile 
(a = 1) which we assumed for analysis in $3] The horizontal bar denotes the la statistical 
error on a derived from the probability distribution. 



